Phase transitions induced by saddle points of vanishing curvature 
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Based on the study of saddle points of the potential energy landscapes of generic classical many- 
particle systems, we present a necessary criterion for the occurrence of a thermodynamic phase 
transition. Remarkably, this criterion imposes conditions on microscopic properties, namely curva- 
tures at the saddle points of the potential, and links them to the macroscopic phenomenon of a phase 
transition. We apply our result to two exactly solvable models, corroborating that the criterion de- 
rived is not only valid, but also sharp and useful: For both models studied, the criterion excludes 
the occurrence of a phase transition for all values of the potential energy but the transition energy. 
This result adds a geometrical ingredient to an established topological condition for the occurrence 
of a phase transition, thereby providing an answer to the long standing question of which topology 
changes in configuration space can induce a phase transition. 

PACS numbers: 05.70.Fh, 05.20.-y, 75.10.Hk 
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A phase transition is an abrupt change of the macro- 
scopic properties of a many-particle system under varia- 
tion of a control parameter. Typical examples of phase 
transitions are the sudden disappearance of the electric 
resistance when cooling a superconducting material be- 
low its transition temperature, or the evaporation of a 
liquid at temperatures above its boiling point. Consid- 
erable progress has been made in the understanding of 
various aspects of phase transitions, in particular in the 
study of critical phenomena by renormalization group 
techniques. Other aspects are less understood. One of 
those is the question under which conditions on, say, the 
Hamiltonian of a system we can expect a phase transi- 
tion to occur, and when can we exclude such a transi- 
tion? Very little is known about conditions guaranteeing 
a phase transition to take place. The Peierls argument [1| 
or the Frohlich-Simon-Spencer bound [H can be used to 
prove the existence of phase transitions without explicitly 
computing a thermodynamic potential, but the applica- 
tion of these methods is model specific and, depending 
on the system of interest, may be very difficult. There 
are a few remarkable results giving necessary conditions 
for the occurrence of a phase transition, among those 
the Mermin- Wagner theorem 0] or the theorems on the 
absence of phase transitions in certain one-dimensional 
systems by van Hove Q and by Cuesta and Sanchez 
Such results are very useful since they can exclude the oc- 
currence of phase transitions not only for specific models, 
but for classes of systems. 



One attempt to better understand the origin of phase 
transitions and the conditions under which they can oc- 
cur dates back to Caiani et al. @. These authors con- 
jectured that phase transitions in classical many-particle 



systems were related to topology changes of the surfaces 



£„ = {q G Tjy | V(q) = Nv} 



(1) 



of constant potential energy (per particle) v = V/N in 
configuration space T jv , where N is the number of degrees 
of freedom. In the ensuing years, this "topological hy- 
pothesis" was corroborated by a number of model calcu- 
lations, although some peculiar (e.g., long-range) models 
disprove the all-encompassing validity of the conjecture 
(see [7[ for a review). Only recently, a connection between 
phase transitions and topology changes in configuration 
space was demonstrated rigorously: For a certain class of 
systems with short-range interactions, a topology change 
of at v — vt was shown to be necessary for a phase 
transition to take place at a transition potential energy vt 
This theorem, at least for the systems satisfying its 
assumptions, suggests the interpretation of certain topol- 
ogy changes as the relevant mechanism for the occurrence 
of a phase transition. Furthermore, in the manner of the 
other necessary conditions mentioned above, the theorem 
allows one to exclude a phase transition whenever such 
topology changes are absent. 

Unfortunately, in many realistic situations an applica- 
tion of the theorem in H E| is impossible: For the typi- 
cal systems of interest, the number of topology changes 
is very large, growing unboundedly with N, and in the 
thermodynamic limit N — > oo the topology changes be- 
come dense on a part or even the entire energy axis. So 
most of the topology changes turned out to be unrelated 
to phase transitions, but a criterion to identify the rele- 
vant ones was lacking. In this Letter, we supplement the 
topological approach to phase transitions with a geomet- 
rical ingredient which allows us to quantify the contribu- 
tion of each topology change to thermodynamic quanti- 
ties. Based on a finite-system result put forward in [icj |. 
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we derive a topological-geometrical criterion by which we 
identify most of the topology changes as being unrelated 
to phase transitions in the thermodynamic limit. We 
apply this criterion to two exactly solvable models, the 
mean-field XY model and the mean-field fc-trigonomctric 
model. For both models studied, our criterion excludes 
the occurrence of a phase transition for all but one value 
of the potential energy, coinciding precisely with the po- 
tential energy at which a phase transition occurs. 

Under rather mild conditions on the many-particle po- 
tential V, topology changes of the subsets E„ are closely 
related to saddle points of the V. This connection links 
the topological approach to phase transitions with the 
study of energy landscapes, an active field of research 
with many applications to clusters, bio-molecules, and 
glass- formers [ll]. We use the language of energy land- 
scapes in the following exposition. 

Preliminaries. — We consider classical systems of N de- 
grees of freedom, characterized by a Hamiltonian func- 
tion of standard form, 

1 N 

H(p,q) = -Y,P 2 i+V(q), (2) 

i=l 

where p — (p%, . . . ,pn) and q — (q%, . . . , gjv) are the vec- 
tors of momentum and position coordinates. The restric- 
tion to such a standard form is not essential, but simpli- 
fies the discussion: Since a quadratic form in the mo- 
menta pi yields only an irrelevant additive contribution 
to the free energy, it is sufficient to focus on the potential 
V in the following. The potential V is an analytic map- 
ping from the continuous configuration space T pj onto the 
reals, and in general V has a number of critical points (or 
saddle points) q c with vanishing differential, dV(q c ) = 0. 
We assume in the following that all critical points of V 
have a non-singular Hessian S)v- In this case, V is called 
a Morse function, and one can argue that the restric- 
tion to this class of functions is not a serious limitation: 



Morse functions form a dense, open subset of all analytic 
functions, and hence are generic. Therefore, even if V is 
not a Morse function, it can be made into one by adding 
an arbitrarily small perturbation. A central proposition 
of Morse theory states that, for potentials V being Morse 
functions, each topology change of the subsets 

M v = {q€T N \ V(q) < Nv} (3) 

at some value v = v c corresponds to one or several critical 
points q c of V with critical value v c — V{q c )/N . The type 
of the topology change is determined by the index of a 
critical point, i.e., the number of negative eigenvalues 
of the Hessian Sjy at q c . As a consequence, the study 
of critical points and critical indices of V is equivalent 
to the study of topology changes of M. v . Since E„ as 
defined in ^ is the boundary of Ai v , both subsets are 
closely related. We now continue the exposition in the 
language of saddle points and reinterpret our results in 
the language of topology changes only in the summary. 

Finite system nonanalyticities. — In a recent Letter 
[Io| . we presented an exact, quantitative account of the 
contribution of a single critical point of a Morse function 
V to statistical mechanical quantities of finite systems 
like the (configurational) density of states 

(l N (v) = f dqS[V(q)-Nv] (4) 
Jr N 

or, equivalently, the entropy sn(v) = ln[fijy(w)]/iV per 
degree of freedom. We found that the density of states 
51 at = + fi^f can be written as the sum of a "harm- 
less" analytic term S7jy and a nonanalytic term which, to 
leading order, is given by 

^» = ( N ^) N/2 =h n N \(v) (5) 

with the universal function 



r (_l)fe/2 v (N-2)/2 for k evenj 

h U N, k {v) = \ (-l)( fe +l)/ 2 v (N-2)/2 n -l m | w | for N k oddj (g) 

(_ 1 )(JV-fc)/2(_ u )(A'-2)/2 Q ^_ w j for ^ k od(L 

I 



Here k denotes the index of the critical point and is the 
Heaviside step function. This result, for its relation to the 
index of the critical point, is topological in nature, but it 
has also a geometrical ingredient: the determinant of the 
Hessian in the denominator of ([5|) can be interpreted as 
some sort of curvature, quantifying the "flatness" of the 
critical point. 



Thermodynamic limit. — We now take the above finite- 
system result as a starting point to compute the contribu- 
tion of critical points of V to thermodynamic quantities. 
To this purpose, we define the e-entropy in the thermo- 
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dynamic limit as 

s^(v) = lim 1 In [A^OO + B#'»] (7) 
= max{a Vo ' e (v),b v °' e (v)}, (8) 

where i?^' e contains the nonanalytic contributions JTjy gc 
from the critical points g c in the e- neighborhood of vq , 

BT e (v)= ]T E ^«c(«). ( 9 ) 

{u o :|l> o -u |<e} {<?c:V(g c )=Afu c } 

for some small e > 0. The sum is written as a sum over 
critical values v c of V and one over critical points q c {v c ) 
of the respective critical value. Then a smooth function 
A V j^' e can be chosen such that the e-entropy s Vo,e coin- 
cides with the exact entropy s = limjv-^oo N~ x lnfijv in 

I 

^(^i^{/:>«'<-') ( "->o(»') 



the interval (vo — e,vo + e). From the results in [j| 0| 
one can argue that, at least for well-behaved short-range 
systems, the analytic contributions in A v ^' e are irrelevant 
for phase transitions, converging uniformly to an analytic 
function a Vo,e = limjv^oo ^V" 1 In A^' e in the thermody- 
namic limit. Hence, we examine under which conditions 

b V0 ' e (v) = lim ±hxB%> e (v) (10) 

N~+oo IX 

may contribute to the entropy in the thermodynamic 
limit, thereby possibly inducing a phase transition. 

Inserting the functional form of {Yff into (J9j> , we have 
to distinguish the various cases in Eq. ([6]). For N = 1 
(mod 4), we obtain 



N 2 {v')] + / dv' (v' - v ) {N - 2)/2 [Mi(v') -M 3 (v')] , 

J v J 

(in 



but other cases look similar and yield identical results in 
the thermodynamic limit. Here, 

Nt{v') = J (<lc) S[V(q c )/N - v'] S iMqc) mod 4 (12) 



are the distribution functions of the critical points q c with 
index k(q c ) = £ (mod 4), weighted by their Jacobian de- 
terminant 



J(g c ) = |det[io y (< Zc )/2]r 1/2 . 



(13) 



Under suitable assumptions on the thermodynamic limit 
properties of the relevant quantities we can evaluate the 
integrals in (fTTj) by Laplace's method and construct the 
upper bound 

1 



b v °> e (v) < - m (e) + V27re+ max [m + Mv')] , (14) 
2 e, |u--u'|<£ 

where exp(A r n^) is the overall number of critical points 
with index k = £ (mod 4), generically growing exponen- 
tially with the number N of degrees of freedom. The cor- 
responding Jacobian density in the thermodynamic limit 
is given by 



qc£Qe(v,v+e) 



9c &Qe (v,v+e) 



(15) 

where Qg(v,v + e) denotes the set of critical points q c 
with index k(q c ) = £ (mod 4) and with critical values 
V{q c )/N in the interval [v, v + e]. 



Now, if the max-term in (|14p is finite, we can always 
choose e sufficiently small such that the bound for b v °' e 



is smaller than 



Then, as a consequence of ©, the 



entropy s in the interval (vq — e, vq + e) is exclusively 
determined by a" 0,e , which we assume to be a smooth 
function. This argument can be repeated for any value 
of vq, thereby demonstrating that the saddle point con- 
tributions in the 6-term are harmless in the sense that 
they do not lead to nonanalyticities in the entropy. This 
reasoning, however, relies on the assumption that the 
max-term in (1141) is finite, and it breaks down in case ni 
or je diverge. These observations lead us to the central 
result of this Letter. 

Theorem. The saddle point contribution b Vo,e (v) cannot 
induce a phase transition at any potential energy in the 
interval (vq — e, vo + e) if 

1. the number of critical points is bounded by 
exp(CiV) for some C > and 

2. the Jacobian densities have a thermodynamic limit 
of the form [Ijty with j e < oo W € {0, 1, 2, 3} 

inside the given interval. 



In the literature on energy landscapes [ll|, the first 
condition is assumed to be true for generic many-particle 
systems, although an actual proof is missing. The second 
condition in the Theorem is the key result which allows 
us to identify most of the saddle points of the potential 
V, and therefore most of the topology changes in config- 
uration space, as irrelevant for the occurrence of a phase 
transitions. Only the ones which are sufficiently "flat" 
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such that the Jacobian density je diverges in the ther- 
modynamic limit can give rise to a phase transition. To 
corroborate our reasoning and to illustrate the content of 
the Theorem, we present the computation of je for two 
simple, exactly solvable spin models, showing that phase 
transitions occur precisely at the values of v for which 
condition [2] of the Theorem breaks down. 

Examples. — The mean-field XY model is a system of 
N plane rotators described by angular variables qi G 
[0, 2w) (i = 1, • • • , N). Each rotator is coupled to each 
other with equal strength K > 0, where the interactions 
are described by the potential 



V XY (q) 



K 
2N 



N 



N 



[1 - cos(^ - gj)] - h' s ^ j cos q L . (16) 



i=l 



The rotators are subject to a magnetic field of strength 
h which energetically favors orientations qi ~ it for h < 
and qt ss (mod 2ir) for ft, > 0. In the limit N — > oo and 
h — > 0, the system shows a continuous phase transition. 

The critical points of this model and their indices have 
been analyzed by Casetti et al. [l2| for arbitrary numbers 
N of rotators. Critical points of Vxy{q) were found to 
occur for q = q c E {0, 71-}^, i. e., for (q l5 . . . , qjy) with all 
components qi being either or n. Hence, in accordance 
with condition [T] of the Theorem, the number of critical 
points is growing exponentially with N. The elements 
pa = d 2 V/{dqidqj) of the Hessian are also reported in 
[12j| , Evaluated at a critical point, the diagonal elements 
of the Hessian can be written as 



K 1 



COS qi 



K 



(17) 



where is the number of 7r's in the sequence 
(<7ij ■ • ■ i On) & {0,tt} n ■ It has been shown in Appendix 
A.l of [12] that, in the thermodynamic limit N — > oo, the 
contribution of the off-diagonal elements to the eigenval- 
ues of 5) vanishes. Therefore, for large N we can write 
the Jacobian determinant (|13|> at a critical point q c as 



jv»i 



N 



Y[Sju(q c ) 



i=l 



-1/2 



K 
~2 



N 



-N/2 



(18) 



Inserting an expression of depending on the potential 
energy per degree of freedom v, 



N 



1 -*jx( K + h± ^ R2 + h2 - 2Kv ) . ( 19 ) 



we can calculate the Jacobian density ([15| for N — > oo, 

£-0,1,2,3, (20) 



jt(v) = ±ln2-^ln[K(K-2v) 



in the limit h — * 0. The graph of this function is plotted 
in Fig. [1] (left). As is easily recognized, je(v) is singular 
only at v = K/2, which is precisely the potential energy 



per degree of freedom at which a phase transition occurs 
for the mean-field XY model with zero external field. In 
agreement with our Theorem, no phase transition is in- 
duced by the critical points with finite Jacobian density. 

The mean-field fc-trigonometric model, characterized 
by the Hamiltonian function 

N 

V k (q) = AN 1 - k Yl [1 + "• + »»)]. (21) 



is another spin models consisting of plane rotators, with 
A determining the coupling strength between the rota- 
tors. For this model, the thermodynamic behavior and 
the critical points of the potential have been studied by 
Angelani et al. [13| . The parameter k € N in the Hamil- 
tonian function crucially determines the thermodynamic 
behavior of the model: In the limit N —> oo, the system 
shows a discontinuous phase transition for k 3, a con- 
tinuous one for k = 2, and no phase transition for k = 1. 
The results on the critical points of Vk in [13| can be as- 
sembled, in a calculation very similar to the one for the 
mean-field XY model, to obtain the Jacobian densities 
ji as functions of the potential energy per particle v, 



i/fe-i 



£=0,1,2,3. (22) 



The graph of this function is plotted in Fig. [T] (right) 
for k — 1,2,3,4. je is constant (and therefore bounded 
above) for k = 1, in agreement with the absence of a 
phase transition in this case. For k 2 the Jacobian 
density ji(v) shows a divergence at v = A, coinciding 
precisely with the transition potential energy of the phase 
transition of the mean-field fc-trigonometric model and il- 
lustrating nicely the content of the Theorem. In contrast 
to short-range systems, phase transitions in long-range 
systems need not be, but can be induced by saddle points 
of V with diverging jj>. According to our results, this is 
the relevant mechanism for the two models studied. 

Summary. — Based on a result describing the effect of 
saddle points of the potential V on the finite-system en- 
tropy, we derived a necessary criterion for the occurrence 



A 



FIG. 1: (Color online) Plot of the graphs of the Jacobian 
densities je as functions of the potential energy v for the 
mean-field XY model with coupling K = 1 (left) and for 
the mean-field fc-trigonometric model with k = 1,2,3,4 and 
A = 1 (right). In agreement with our Theorem, the transi- 
tion potential energies vt ~ 1/2 and v t — A of these models 
coincide with the singularities of je. 
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of a phase transition in the thermodynamic limit. This 
criterion imposes conditions on microscopic properties, 
namely curvatures at the saddle points, and links them 
to the macroscopic phenomenon of a phase transition. 
Interpreting the result (via Morse theory) in terms of 
topology changes in configuration space, an answer is 
provided to the long standing question of which topol- 
ogy changes of or M. v can induce a phase transition. 
Applications to the mean-field XY and fc-trigonometric 
models show that the criterion is not only valid, but also 
sharp and useful: For both models, the criterion excludes 
the occurrence of a phase transition for all values of the 
potential energy but the transition energy. 

Our Theorem provides a criterion to single out those 
saddles of the potential energy which are relevant for 
phase transitions, and it facilitates the study of phase 
transitions from the alternative point of view of energy 
landscapes. This approach is of particular interest for 
applications to the folding transition of proteins: energy 
landscape methods are frequently used for studying dy- 
namical properties of such transitions, and our results 
allow investigation of both, dynamical and thermody- 
namical properties within a unified framework based on 
energy landscapes. 
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